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Abstract 

> 

0^ . It is well established that the onset of spatially periodic vortex states in the Taylor- 

. Couette flow between rotating cylinders occurs at the value of Reynold's number 

. predicted by local bifurcation theory. However, the symmetry breaking induced by 

ff^ ■ the top and bottom plates means that the true situation should be a disconnected 

■ pitchfork. Indeed, experiments have shown that the fold of the disconnected branch 
can occur at more than double the Reynold's number of onset. This leads to an 

Ch ■ apparent contradiction: why should Taylor vortices set in so sharply at the value 

^ ■ Reynold's number predicted by the symmetric theory, given such large symmetry- 

^ ' breaking effects caused by the boundary conditions? This paper offers a generic 

■ explanation. The details are worked out using a Swift -Hohenberg pattern formation 
model that shares the same qualitative features as the Taylor-Couette flow. Onset 
occurs via a wall mode whose exponential tail penetrates further into the bulk of the 
domain as the driving parameter increases. In a large domain of length L, we show 
that the wall mode creates significant amplitude in the centre at parameter values 
that are 0(L~^) away from the value of onset in the problem with ideal boundary 
conditions. We explain this as being due to a Hamiltonian Hopf bifurcation in 
space, which occurs at the same parameter value as the pitchfork bifurcation of the 
temporal dynamics. The disconnected anomalous branch remains 0(1) away from 
the onset parameter since it does not arise as a bifurcation from the wall mode. 
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1 Introduction 

The Taylor-Couette experiment provided one of the first quantitative verifi- 
cations of the correctness of the Navier-Stokes partial differential equations 
(PDEs) describing the dynamics of fluid flows. The experiment, in its simplest 
form, consists of a pair of concentric cylinders with a fluid-filled gap in be- 
tween; as the inner cylinder is rotated, a shearing flow (the Couette flow) is 
established between the cylinders, and this becomes unstable to axisymmetric 
vortices (Taylor vortices) at a critical value of the rotation rate (as measured 
by a dimensionless Reynold's number R). See [1,2] for reviews. One notable 
achievement of Taylor's work [3] in 1923 was the theoretical prediction and 
the experimental measurement of the critical Reynold's number Rc for the 
onset of vortices, with remarkably good agreement between the two. In doing 
the stability calculation, Taylor assumed that the vortices would be periodic 
in the direction along the axis, and neglected the effects of the top and the 
bottom plates of the experiment. With this assumption, the bifurcation lead- 
ing to Taylor vortices is a pitchfork (figure la), and the characteristic sharp 
transition of this bifurcation, with the strength of the vortices going as the 
square root of the degree of supercriticality R — Rc, has been confirmed exper- 
imentally [4]. The sjTiimetry that is broken in this pitchfork is a translation 
symmetry along the axes of the cylinders. 

Subsequent theoretical developments explored the role of the top and bottom 
plates in the experiment, which spoil the idealisation of spatial periodicity in 
the direction parallel to the axis of rotation, and which break the translation 
symmetry assumed in the original theoretical work. Ekman boundary layers 
cause the fluid near the boundary to spiral preferentially inwards for any non- 
zero rotation rate, and the Taylor vortices to develop first in the boundary 
layer, moving smoothly into the bulk of the fluid as R approaches Rc - this 
has been observed in experiments [7] and in calculations [8]. In the words of 
Benjamin [5], 'no precise critical value of R exists for the onset of cellular 
motion'. Benjamin [5] interpreted the formation of Taylor vortices as a broken 
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Fig. 1. Sketches of the amphtude of Taylor vortices measured in the centre of the 
apparatus as a function of Reynold's number i?, under various assumptions, (a) With 
ideal (reflecting or periodic) boundary conditions, there is a sharp transition to 
Taylor vortices at a pitchfork bifurcation &t R = Rc- (b) With the assumption 
of a weakly broken pitchfork bifurcation, there is still a relatively sharp onset of 
Taylor vortices close to R = Rc (after [5]); (c) The experiments of Benjamin and 
Mullin [6] suggest that the upper half is still a weakly broken pitchfork, whereas 
the lower half is a strongly broken pitchfork, with anomalous modes only appearing 
at Reynold's numbers at least twice the critical value. Solid (dashed) lines indicate 
stable (unstable) solutions. 
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Fig. 2. Finite-element results reproduced with permission from [9] showing 
two radial velocity profiles below and above the critical value of R at which 
the onset of periodic vortices would occur with periodic boundary conditions, 
(a) R/Rc = 63.74/68.189 = 0.935; (b) R/Rc = 69.02/68.189 = 1.013. 



pitchfork (figure lb), with the end plates driving a flow near the boundary 
for all non-zero R, and this flow exciting a cellular flow that penetrates the 
central region with increasing Reynold's number - see figure 2. 
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The qualitative role of the end boundaries was explored further by Schaef- 
fer [10], who introduced a homotopy parameter r (0 < r < 1), with r = 
corresponding to ideal end boundary conditions, where a state of pure Cou- 
ette flow exists for all values oi R, t = 1 corresponding to physically realistic 
end boundary conditions, and intermediate values of r interpolating between 
these two extremes. The results for r close to zero are indeed consistent with 
a weakly broken pitchfork bifurcation (figure lb), and are in qualitative agree- 
ment with experimental results of Benjamin [5,11]: as the Reynold's number 
is increased slowly, vortices grow smoothly, with the most rapid growth occur- 
ring for R close to the critical value. This is illustrated in figure 2: note how 
the radial velocity profiles are not zero for R below Rc] instead we see a pair 
of 'wall modes' connecting the non-parallel flow at the walls to Couette flow 
in the main body of the cylinder. Similarly, the steady solution for R greater 
than Rc is not a pure periodic solution but has modulation near the two walls 
in order to satisfy the inhomogeneous boundary conditions. 

However, by starting the experiment impulsively, it is possible to find another 
branch of vortices {anomalous modes) that resemble ordinary vortices, but 
have the opposite sign - close to the ends, these anomalous vortices have an 
outwards radial velocity, opposite to the normal vortices [11]. Other experi- 
ments and calculations indicate that anomalous modes may also have stag- 
nation regions or narrow counter-cells close to the top and bottom bound- 
aries [12,13]. 

However, anomalous modes cannot be found close to the critical Reynold's 
number, as would be suggested by figure lb. Instead, R needs to be at least 
twice its critical value (and the cylinder started impulsively) before anoma- 
lous modes can be found [6] . Once they are established, the anomalous modes 
persist as R is decreased to a lower stability bound, the exact value of which 
depend on experimental parameters such as the gap width between the in- 
ner and outer cylinders, or the aspect ratio L, a dimensionless measure of 
the length of the column. The lower existence boundary appears always to 
be at least twice the critical value Rc [6], and can be much higher if the gap 
between the cylinders is narrow [14]. Thus the experimental situation is de- 
picted in figure Ic. Interestingly, the lower stability boundary of the anomalous 
modes seems to be independent of the aspect ratio of the apparatus (for large 
aspect ratios), and remains at an appreciable multiple of the critical Reynold's 
number for ordinary modes even as L becomes large [14,15]. This emphasises 
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the fact that the large (but finite) aspect ratio hmit is very different from the 
ideahsation of periodicity in the axial direction. 

It is worth emphasising the two surprising and apparently contradictory as- 
pects that have emerged. First, the distance {d in figure 1(c)) between the fold 
on the disconnected (anomalous) branch is such that the Reynold's number 
at the fold is at least a factor of two larger than the Reynold's number of 
the onset of ordinary Taylor vortices. Making the cylinder longer (so that the 
boundary effects are moved 'towards infinity') does not make d tend to zero, 
and many authors have concluded that the onset of Taylor vortices can in no 
way be regarded as a weakly broken pitchfork. Second, the onset of vortices 
is sharp when viewed in terms of measures such as the radial velocity of the 
midpoint of the apparatus, and it occurs at almost exactly the value of R 
that is predicted for the problem without end effects - so the onset of Taylor 
vortices apparently can be described as a weakly broken pitchfork bifurcation. 
This discrepancy would not be expected for a generic unfolding of a pitch- 
fork bifurcation. The purpose of the present paper is to explain this apparent 
contradiction. 

The rest of the paper is outlined as follows. In section 2, we introduce the 
Swift-Hohenberg equations as a model for the Taylor-Couette flow. We note 
that the modiflcation of pattern formation due to the presence of weak forc- 
ing at lateral boundaries in Swift-Hohenberg equations has been addressed 
in the work of Daniels and co-workers [16,17] with application to Rayleigh- 
Benard convection in mind. In contrast with their work, we include an order 
one boundary condition that forces the flow strongly. This approach is also 
distinguished from that of [18], which focused on a Ginzb erg-Landau equation 
for the envelope of the vortex amplitude. 

Section 3 then includes an analysis of the linearised version of this model. 
Despite being purely linear, it is found that mode shapes like those in flgure 2 
emerge under static increase of the bifurcation parameter through the value 
at which the symmetric problem bifurcates. We explain this in terms of the 
relation between the temporal pitchfork bifurcation and a spatial Hamiltonian 
Hopf bifurcation. 

Section 4 goes on to consider a nonlinear bifurcation analysis. It transpires 
that the unfolded pitchfork resembles flgure 1(c). There are many anomalous 
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branches that emerge from primary and secondary symmetry breaking bifur- 
cations in the symmetric problem. In the Swift-Hohenberg example studied 
in detail, the stable anomalous branch does not always emerge from the pri- 
mary pitchfork bifurcation when the symmetry-breaking terms tend to zero. 
We explore this issue in some detail. Finally, section 5 draws conclusions and 
discusses wider implications of the results. 

2 Swift— Hohenberg model 

Rather than consider the nonlinear axisymmetric hydrodynamic partial differ- 
ential equations (PDEs) that describe the flow between two rotating cylinders, 
we focus on the simpler Swift-Hohenberg [19] PDE, which shares many of 
the same pattern- forming features. In fact, Melbourne [20] has demonstrated 
that bifurcation problems of the Taylor-Couette type (steady state bifurca- 
tions with nonzero critical wavenumber in systems with Euclidean symmetry) 
reduce to equations of Swift-Hohenberg form (though with more general non- 
linear terms). The model equation we use is: 

Ut = fiU-il + 2U,,^ + U^^^^) -U^-U U^, (1) 

where the subscripts denote partial derivatives with respect to space x and 
time t. The dependent variable U{x,t) G M is defined on x G [— L/2,L/2], 
where L represents the length of the Taylor column. The parameter /i repre- 
sents the forcing R — R^. The usual form of the Swift-Hohenberg equation has 
only a cubic nonlinearity, but we include a quadratic term to ensure there is no 
U —U symmetry (see, e.g., [2, eq. (7.21)]. This term also makes the model 
equation non-variational and so allows unsteady behaviour as an asymptotic 
state (though we focus entirely on steady states). 

In order to relate the order parameter U to the fluid flow, we interpret f/ as a 
stream function, and so represents the radial velocity in the column. The 
effect of the non-slip boundary conditions on the top and bottom plates is to 
induce an inwards flow near the boundaries, although the radial velocity is 
zero on the boundaries themselves. We model this strong forcing at the end 
walls with the inhomogeneous boundary conditions: 

t/(-L/2) = f/(L/2) = 0, U^{-L/2) = U^{L/2) = -1. (2) 
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The only symmetry in the problem is the reflection in the equatorial mid-plane 
of the apparatus: 

{U,x)^i-U,-x). (3) 



We also consider idealised reflecting boundary conditions: 



U{-L/2) = U{L/2) = 0, 



U,,{-L/2) = U,,{L/2) = 0, 



(4) 



which have an additional hidden symmetry [21]: the problem can be extended 
by reflection onto the domain x G [—L, L] with periodic boundary conditions 
and so acquires a continuous translation symmetry. It is this translation sym- 
metry that is broken in the pitchfork bifurcation in the idealised version of this 
problem (figure la). This symmetry is strongly broken by the inhomogeneous 
boundary conditions (2). 

Experiments find states where the mid-plane reflection symmetry is preserved, 
so we focus on this case by using reflecting boundary conditions at x = 0. 
We are also interested in examining the transition from idealised to realistic 
boundary conditions, and so we use: 



U{0) = U,,{0) = U{L/2) = r {U,{L/2) + 1) + (1 - r)U,,{L/2) = 0. (5) 



Here r(0<r<l)isa homotopy parameter [10], such that r = corresponds 
to the ideal problem (with hidden translation symmetry) and r = 1 to the 
realistic non-slip boundary conditions. 

3 Linear analysis: pitchfork and Hamiltonian Hopf 

Solutions of the linearised problem can be written in terms of exponentials in 
time and space, in the form: 



where s is the temporal growth rate, /c is a spatial wavenumber, and a is a 
spatial growth rate. In order to satisfy the linearised model equation, s, a and 
k must satisfy 



U{x,t) = e 



st+(o'+ifc)a; 



(6) 



= ka{k^ 



2 

— (7 — 



1), 



s = \i — [k 



2 



2 

— O — 



l)2 + 4aV, 



(7) 
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which can be rearranged to give three possibihties: 



B 



C 



A 




s = fi- + 

s = fx - {1- k^, 

s = n + 4o-^(l + o-^) 



(9) 
(10) 
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Since the hnearisation of the PDE (1) is first order in time and fourth order 
in space, there is a unique temporal growth rate s and a total of four complex 
spatial growth rates, which are roots of the equations above, namely in case A: 
±0"! and ±(72; in case B: ±iki and ±ik2] in case C: ±cr ± ik. Corresponding 
to these, the linear solutions are of the form: 



where A and B are constants that will be determined by the boundary con- 
ditions at X = L/2. The odd boundary conditions at x = have already been 
enforced by the choice of trigonometric functions. 

The next stage of the calculation depends on whether ideal (r = 0) or realistic 
(r = 1) boundary conditions are being used. With ideal boundary conditions, 
t/ = is always a solution of (1), and bifurcations from this state occur 
when there are marginally stable (s = 0) linear solutions. Setting s = 
and U{L/2) = Uxx{L/2) = results in no solution in cases A and C, and an 
eigenvalue problem in case B, where k can take on discrete values: k = 27m/ L, 
where n is the number of vortices in the half-domain. There are thus pitchfork 
bifurcations at 



where n is an integer. In particular, there is a pitchfork bifurcation at = 
whenever the domain is chosen to fit an exact number of vortices: n = L/2tt. 
Note that the condition for the onset of the vortices is s = and a = 0, 
with a purely imaginary spatial wavenumber ±ik corresponding to a spatially 
periodic pattern. On either side of this bifurcation point, the temporal growth 
rate s, indicating the stability of the Couette flow, changes from negative to 
positive at the pitchfork bifurcation. 

On the other hand, with inhomogeneous boundary conditions (r = 1), there 



A: U{x,t) = e^* {A sinh(o-ix) + B sinh(cT2x)) , 

B : U{x,t) = e** {A sin(A;ix) + B sm{k2x)) , 

C : U{x,t) = e** {A cos{kx) sinh(cTa;) + B sm{kx) cosh(cTa;)) , 



(11) 
(12) 

(13) 




(14) 
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Fig. 3. Dependence of the radial velocity at the midpoint {Ux{0)) as a function of 
bifurcation parameter fi for the nonlinear (solid) and linearised (dashed) PDE (1), 
with L = 227r and realistic boundary conditions (r = 1). Also shown in (b) as a 
dotted line is an integrated average of the solution and its first three derivatives, 
showing how other measures of amplitude are not so sharp as the radial velocity at 
the midpoint. 



are nonzero steady (s = 0) linear solutions for all values of fi. Case A can 
only arise when fi > 1 (and in fact only a single solution is possible); in 
case B, there are two possible solutions for > /i > 1 and one for /i > 1, 
and case C is possible only when /i < 0. The solutions are thus, for /i < 0: 
U{x) = Acos{kx) sinh((Tx) + B sin{kx) cosh.{ax), where a > and k > are 
determined from fi by 4(t'^ + 4o"^ + /i = and /c^ = 1 + o"^; and for < /i < 1: 
U{x) = 74sin(A;ix) + i?sin(A;2x), where fci > 0, A;2 > and kf = 1 + ^/JI, 
k1 = 1 — ^JJi. In these two expressions for U{x), the constants A and B 
are determined by a pair of linear equations from the boundary conditions 
U{L/2) = and t4(L/2) = —1. The dependence of the solution, as measured 
by the radial velocity at the midpoint {Ux{0) = Aa + Bk for /i < 0, and 
Ux{0) = Aki + Bk2 for < /i < 1) is shown in figure 3 (dashed line). 

The amplitude of the linear solution goes to infinity for /i about 0.00826 
(with L = 227r). Since there is no sharp onset with realistic boundaries, one 
cannot define a precise value of fi at which pattern will be first observed in a 
domain of finite length, but the value of = fioo for which the linear solution 
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goes to infinity is a suitable proxy. This is defined implicitly by the condition 



where ki^2 are defined in terms of /i above. When L is large, the smallest 
positive solution fj,oo occurs when {ki — k2)L/2 = n. This yields 



In other words, the divergence of the linear solution occurs for /i closer to 
zero as L increases. Moreover, we have that for large L the first pitchfork 
bifurcation occurs at 



same order: L . 

Also shown in figure 3 are nonlinear solutions of (1). The linear and nonlinear 
solutions are close to each other for /i negative: the discrepancy at low ampli- 
tude arises because the linearised solution is not small close to the boundaries. 

Linear solutions of the PDE (1) are shown in figure 4 for two values of /i. As /i 
is increased from negative values to = 0, the exponentially decaying linear 
solution extends further into the bulk of the fluid. In an arbitrarily long cylin- 
der, the radial velocity at the midpoint of the apparatus would remain almost 
zero until the spatial decay rate (as measured by a) became zero. Therefore, 
the condition for onset of steady vortices, as measured at the midpoint of the 
apparatus, is a = and s = - which is the same condition as for the onset 
of vortices with the idealised boundary conditions. In this case, however, on 
either side of onset, it is the spatial eigenvalues ±0" ± ik and iz/ci, ±zA;2 that 
change in nature, at a Hamiltonian Hopf bifurcation [22] (see figure 5). 

More accurately, we should describe the spatial bifurcation as a reversible 
1:1 resonance, since the ODE obtained by setting = in (1) does not 
conserve a first integral and so cannot correspond to a Hamiltonian system, 
but is nevertheless reversible in the sense analysed by looss and Peroueme [23] 
owing to the symmetry (3). Here the spatial bifurcation is supercritical and for 
/i > there exist spatially periodic solutions the maximum amplitude of which 
grows as the square root of for the nonlinear problem. The implications of 
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Fig. 4. Linear solutions of (1), with L = 40 and r = 1, with (a) fi = —0.1 and 
(b) n = 0, showing the growth of the solution in from the edges (x = ibL/2). 
Compare with figure 2. 



(a) 



(b) 



(c) 



Fig. 5. Hamiltonian Hopf bifurcation: (a) /U < 0, the spatial eigenvalues are ztcntiA;; 
(b) /X = 0, the eigenvalues ±i; (c) /i > 0, the eigenvalues are ±iki, ±i/c2. 

this bifurcation are discussed in more detail in Section 5. 



This explanation carries over from the Swift-Hohenberg model to the Taylor- 
Couette problem and more general pattern forming situations with strong 
forcing at the boundaries. If one assumes periodic boundary conditions, with 
pure imaginary spatial wavenumbers, the criterion for onset in a general pat- 
tern forming problem is that the temporal growth rate is zero. If, one the other 
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(a) (b) 




Fig. 6. As the length of the domain increases, the pattern, as measured in the centre, 
sets in more sharply: (a) detail near ^ = 0, for L = 227r, 447r and SSvr; (b) the larger 
picture, showing the anomalous modes. 

hand, one takes the end walls into account but the domain is very large, and 
the pattern is measured only far away from the boundaries, then the steady 
inhomogeneous solution will penetrate into the bulk of the fluid and reach the 
centre when the real part of the spatial wavenumber is zero. Thus the two 
perspectives will yield the same condition for the onset of pattern formation: 
yU = in the case of the model PDE, or R = Rc in the case of Taylor-Couette 
flow. This explains the sharp transition seen in large domain Taylor-Couette 
experiments at the Reynold's number predicted using idealised boundary con- 
ditions, even though the boundaries are forcing the flow strongly. 

4 Nonlinear steady-state bifurcation analysis 

The remaining issues to be addressed are the effect of the length of the domain 
on the nonlinear solutions, and the location of the saddle-node bifurcation on 
the anomalous branch (figure Ic). 

We focus on the steady state problem given by setting Ut = in (1). 

Uxxxx + 2Uxx + il-fi)U + U^ + UUx = 0, (18) 
subject to boundary conditions (5). Nonlinear solutions are computed using 
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Fig. 7. As ^ increases through zero, the wall mode penetrates further into the bulk, 
but the amplitude at the centre of the domain remains small for larger values of L. 
(a,b) L = 50, /i = -0.1, 0; (c,d) L = 200, /i = -0.1, 0; (e,f) L = 1000, ^ = -0.1, 0. 
These are nonlinear solutions, but linear solutions look very similar. 

AUTO [24] as a boundary value solver. The effect of the size of the domain is 
illustrated in figures 6 and 7. As the length L increases, the amplitude, as mea- 
sured in the centre, sets in more sharply as /i is increased through zero, and 
the curve resembles half a pitchfork as L ^ oo. With very large values of L 
(figure 7), the exponential decay into the bulk ensures that the inhomogeneous 
pattern has very small amplitude for /i < 0. This is consistent with our under- 
standing from the linear theory. The anomalous mode branches for L = 227r, 
447r and SSvr are also shown in figure 6. Note how the anomalous modes do 
not approach n = for larger values of L, even though the symmetry-breaking 
effects are being pushed further away. 

With a fixed value of L, only one sign of Ux{0) is possible near the transition 
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at /i = 0, but which one is observed depends whether the number of vortices 
between x = and x = L/2 is even or odd. 

Next, we consider the connection between pitchfork bifurcations (in the case 
of ideal boundary conditions) and saddle-node bifurcations (in the case of real- 
istic boundary conditions). With ideal boundary conditions, pitchforks occur 
both from the trivial solution and as secondary bifurcations from the various 
primary branches. We concentrate on the case L = 227r. Figure 8 shows the 
bifurcation diagram computed with ideal (a,b: r = 0) and realistic (c: r = 1) 
boundary conditions. Bifurcation points from the trivial solution (and sec- 
ondary bifurcations from the primary branches) are marked in the figure. The 
primary bifurcation points occur at /i = fi-^f (14), with the bifurcating branch 
being locally proportional to ^fi — /ipf sin(n7r/L). With L = 227r, the solution 
with n = 11 bifurcates precisely from /i = and corresponds to a pattern with 
1 1 vortices in the half domain. The next four bifurcating branches for positive 
fi are also shown in the figure. These bifurcate at /i = 0.0301, 0.0361, 0.1093 
and 0.1574 and correspond to n = 10, 12, 9, and 13 vortices respectively. Note 
that these solutions are invariant under reflections in the midpoint (x = 0), 
and are spatially periodic with period 2L/n. 

Figure 8(c) shows a few of the many branches that exist for realistic boundary 
conditions (r = 1), in the same parameter range as figure 8(b). These were 
obtained by taking all the points on the branches in figure 8(a) with /i = 2 
and continuing these to r = 1, and then continuing in fi once more. Note that 
all of the pitchfork bifurcations have been destroyed, and have been replaced 
by a series of saddle-node bifurcations. The fundamental pitchfork bifurcation 
at /i = has been replaced by a smooth transition, though the remnant of the 
pitchfork shape can clearly be seen, and there is a sharp rise in amplitude close 
to yU = as predicted by the linear theory of the preceding section. Nonlinear 
solution profiles on this fundamental branch (figure 9b) are qualitatively sim- 
ilar to those with ideal boundary condition (figure 9a), and so might be said 
to correspond to = 11 vortices. 

This now brings us to the question of the supposed disconnected part of the un- 
folded pitchfork. Examples of solutions at the labelled saddle-node bifurcation 
points are shown in figure lO(a-h), with the solution Ux{x) at the saddle-node 
bifurcation points drawn as solid lines. These saddle-node bifurcations were 
then continued in (/x, r) back to r = to discover where they originate. Most 
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Fig. 8. (a,b) The idealised bifurcation diagram with r = and L = 22ir: Ux{0) 
as a function of bifurcation parameter /i. Diamonds represent pitchfork bifurcation 
points. Only the first five bifurcating branches from the trivial solution are depicted. 
In increasing order of /i, these have 11, 10, 12, 9 and 13 pairs of vortices in the 
full domain, (c) Bifurcation diagram with r = 1, showing the smooth onset of 
the 11 vortex solution, and several disconnected branches. The plus signs indicate 
solutions that are depicted in figure lO(a-f), at saddle-node bifurcation points. The 
thick lines represent branches that are known to be stable. 
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Fig. 9. Examples of n = 11 vortex solutions with (a) ideal (r = 0) and (b) realistic 
(r = 1) boundary conditions, with fi = 2 (largest amplitude), fi = 1, fi = 0.25, 
fj, = 0.1, /i = 0, ^ = -0.1. 

connect to secondary bifurcations in figure 8(a) and therefore represent un- 
stable solutions (and are shown as dashed lines in figure 10). However, the 
saddle-node bifurcation on the stable anomalous 11 vortex branch, shown in 
figure 10(e) is found to create a stable branch. (The stability of a typical solu- 
tion on this branch was checked by solving, in addition to the ODEs (18), the 
linear variational equations governing by a temporal eigenmode with eigen- 
value s. AUTO was then used to continue solutions in s to establish that there 
are no nontrivial solutions for s > 0.) Unexpected behaviour was found for 
this branch upon varying the homotopy parameter r. One might imagine that 
under homotopy to r = this saddle-node bifurcation point should approach 
the fundamental pitchfork bifurcation at /i = 0. This is not the case for these 
parameter values: when r is reduced from 1 to 0, the saddle-node bifurcation 
itself undergoes a pair of fold bifurcations, and ends up (at r = 0) in the 
unfolding of the n = 12 bifurcation point at /i = 0.0425 (the dashed profile in 
figure 10(e) corresponds to a 12 vortex pattern). 

Alternatively, one could try following the saddle-node bifurcation point that 
occurs as one unfolds the pitchfork at = under infinitesimal increase of r 
from zero. When this is done, the saddle-node bifurcation can be continued 
up to (/i, r) = (4.6220, 1) - see figure 10(h) and figure 11. 

The details of which saddle-node bifurcation (with r = 1) connects to which 



16 



(a) (b) 




X U2 Q X U2 



(g) (h) 




X U2 X L/2 



Fig. 10. Examples of anomalous solutions at the labelled saddle-node bifurcation 
points from figure 8(c). (g,h) are outside the range shown in figure 8. Solid lines 
depict the solution at the saddle node bifurcation with r = 1, and dashed lines 
show the solution that has been continued to r — > 0, ending up at one of the 
bifurcation points in figure 8(a). (a) fi = 0.1367 0.0901; (b) fi = 0.1498 0.0770; 
(c) /i = 0.5193 ^ 0.4452; (d) /i = 0.5581 ^ 0.3787; (e) ^ = 1.6794 ^ 0.0425; 
(f) /i = 1.8227 ^ 0.4452; (g) fi = 2.1318 ^ 0.0770; (h) fi = 4.6220 ^ 0.0000. 
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(a) (b) 




X U2 X U2 



Fig. 11. Continuation of the saddle-node bifurcation: (a) (^,t) = (4.622,1.000), 
(b) (/i,r) = (2.000,0.597), (c) (^,r) = (0.100,0.076), (d) (/i,r) « (0.000,0.000). 

pitchfork bifurcation (with r = 0) was found to depend sensitively on the value 
of L. For example, with L = SStt, the saddle-node bifurcation on the 44 vortex 
anomalous branch does continue down to the primary pitchfork to 44 vortices 
at /i = 0, though in this case the anomalous branch is not stable. The details of 
how the branches connect also depends on the particular choice how the ideal 
and realistic boundary conditions are combined via homotopy. For example 
replacing the final in (5) with a '— ' lead to significantly different results. 



5 Conclusion 

The apparent contradiction described in the introduction is resolved, at least 
in the context of the Swift-Hohenberg model considered here. The onset of 
Taylor vortices is not a weakly broken pitchfork bifurcation, owing to strong 
inhomogeneous boundary forcing. Anomalous modes stay bounded away from 
the critical value of the bifurcation parameter as they must overcome the 
strong preference set by the boundary. When the amplitude of the pattern is 
measured far away from the boundaries, the pattern appears to set in sharply. 
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in half a pitchfork bifurcation, as the decaying wall mode penetrates the bulk 
of the domain. The parameter value at which pattern, as measured in the 
centre of a large domain, become significantly different from zero is the same 
as the value predicted assuming idealised boundary conditions, because the 
requirements for both situations are the same: steady (s = 0) and zero spatial 
growth rate (cr = 0). Moreover we have shown that for a long but finite 
domain of length L, that the parameter value corresponding to this large 
central growth in pattern occurs according to linear theory at a value that is 
within 0(L~^) of the idealised pitchfork. 

The ideas were developed for the model equation, but they apply equally well 
to the Taylor-Couette case, and resolve the difficulties raised by Benjamin 
and Mulhn [6]. 

We have also observed that the saddle-node bifurcation on the anomalous 
branch does not necessarily connect to the primary bifurcation as boundary 
conditions vary from real to ideal. The specific results described here apply 
only to the Swift-Hohenberg model, though the general conclusion that we can 
make is that, under small perturbations from the ideal boundary conditions, 
we expect the pitchfork to be perturbed in the generic way, as in figure 1(b). 
However, going all the way to r = 1 is not a small perturbation, and in general 
nothing can be said about whether the the saddle-node bifurcation created in 
the unfolding of the primary pitchfork is the same (or not the same) as the 
saddle-node at the end of the stable anomalous branch, or even if there is an 
anomalous branch that is stable. 

This approach yields results that are applicable to other pattern formation 
problems (for instance, Rayleigh-Benard convection). Earlier work on con- 
vection [16,17] has focused on weak forcing at the side walls, primarily using 
Swift-Hohenberg theory. We have shown here how the ideas can be extended 
to strong forcing. 

An interesting aspect of our work has been to link the mode selection prob- 
lem to the existence of a spatial Hamiltonian Hopf bifurcation for the infinite 
length problem. This gives the possibility of the existence of branches of spa- 
tially periodic solutions beyond the critical parameter value for the onset of 
rolls, whose period for small amplitude is that given by the wave number of the 
neutral mode of the temporal problem. The normal form of the Hamiltonian 
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Hopf bifurcation (even for reversible systems) is completely integrable up to 
any order (see e.g., [23]) and instead of a unique spatially periodic solution, 
there is a one parameter band of spatially periodic solutions, whose envelope 
grows as the square root of the bifurcation parameter /i. There also exists a 
two-parameter family of spatially quasi-periodic solutions whose existence is 
bounded by the periodic solutions and homoclinic connections to them. How- 
ever, in the reversible case, not all these solutions will necessarily exist in a full 
unfolding of the normal form that breaks its Hamiltonian structure. We note 
from our numerical results for the nonlinear problem with realistic boundary 
conditions and fixed L, that the main branch seems to develop a connection 
from the boundary to a pure periodic state in the middle of the domain. It 
is not clear a priori why this solution and not others are selected from the 
unfolding of the Hamiltonian Hopf bifurcation. 

Finally, we mention another outstanding issue in the problem of the onset of 
Taylor vortices [9]: the timescale for the onset and decay of the pattern are 
different. In particular, these authors computed steady solutions at Reynold's 
numbers just above and just below critical, and in each case altered the 
Reynold's number to an intermediate value and examined the transient. In 
the case of onset, the pattern invaded the bulk as a front travelling in from 
the boundary, while in the case of decay, there was uniform decay throughout 
bulk. The timescales for these two processes were different, and it may be 
possible to explain this using a Swift-Hohenberg based model, as considered 
here. 
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